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Nuclear power plants are often located near a large lake so that cold water can be with- 
drawn to cool the engines. If the lake is thermally stratified so that there is apprciable 
temperature gradient vertically. Cold water can be withdrawn from lake bottom, and warm 
water from the power plant can be returned to the top of lake. Due to the fact that stratifi- 
cation supresses vertical motion, water motion hence thermal mixing should be limited to a 
vertically thin layer not far away from the intake. 

In this section we treat the slow and steady flow of an isothermal but stratified fluid, due 
perhaps to salinity variation, into a two- dimensional line sink. The molecular diffusivity of 
salt is ignored. These simplifications permits a analytical solution which is easy to examine 
the physical implicaions. 

For a slow flow of a viscous fluid flowing in a vertically stratified fluid, we expect, in light 
of Yih's theorem, that the motion of the fluid should be confined within a thin layer. 

2.8.1 Estimation of scales 

Governing equation for continuity : 



u x + w z = 



8.1) 



Ignoring mass diffusion, the incompressiblity condition reads 



up x + wp z = 



(2 



8.2) 



p (uu x + wu + z) 
p (uw x + ww z ) = - 



= -Px + At {u xx + u zz ) 
Py - P9 + A* ( w xx + w zz ) 



(2 
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8.3) 
8.4) 
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Anticipating a thin boundary layer due to vertical suppression, we introduce the following 
normalization. 



then 



US 

u = Uu , w — ~jj~ w ' i P = PoP \ x — Lx , z = 5z 



u' x + w' z = 
u'p' x + w'p' z = 



p U 2 , , , , , Pdp' uU ( , 5 2 

— \p( UUx + WUa )] = - i; — + - 



P U 2 5 . , , , Pdp 1 , U 5 

[p (uw x + w w z )\ = --—- p gp +p^j 



L L 



5 2 L 



U>xxtf+ W zz 



(2.8.5) 

(2.8.6) 
(2.8.7) 

(2.8.8) 
(2.8.9) 



Here 5 is the boundary layer thickness, L ~ x is the horizontal length scale, and U is the 
horizontal velocity scale which must be of the order ~ Q/8. 

After dividing the x,y momentum equations by pU5 2 , and noting 
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PolP f_ = (poUS\ r , 

L pU \ p ) L L 



we get 



Assume 



We must have 



R T \ft Wu' x + w'w' z )\ 



p' x + 



+ u' 



R 



PoUS 



0(1), and - < I. 

Lj 



P5_5_ 
pUL 



= 0(1) 



(2.8.10) 



(2.8.11) 



(2.8.12) 



(2.8.13) 



so that the viscous stress is forced by the pressure gradient. Thus we can take the pressure 
scale to be 

P = P^f (2- 8 - 14 ) 



To the leading order the x momentum equation is simply 

= -p x + pu zz 

Now the z momentum equation, 
PoU 2 5 r , . , , 



L L 



[p (u'w' x + w'w' z )\ 



P 1 dp' , pU 5 

L6/Ldz^~ Po9P+ l 2 ~L 



W xxj2 + W zz 



(2.8.15) 



(2.8.16) 
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Dividing by P/L and using (2.8.14), 



r ( y ] w w + 



i ay p„ 9 l 



(f) 



<9z 



+w ' zz 



i.e., 



I - I [p (w ^ + w wj] = -— - —p + 



Since gravity must be important, we must have 

PogS 



0(1), 



The 2 momentum equation reduces to 



<9p 

= -TT- P9 

oz 



meaning that pressure is hydrostatic. 
Now (2.8.19) implies 



since Q = 0(U5). It follows that 



i.e., 



pogf_ 

pXJL 

pog$ 4 

p,Q L 



0(1) 



0(1) 



4 p,QL 
o ~ 

pog 



Since L ~ x, we have 



' p,Qx 

. pog , 



1/4 



(2.8.17) 



(2.8.18) 



(2.8.19) 



(2.8.20) 



(2.8.21) 

(2.8.22) 
(2.8.23) 

(2.8.24) 



Thus by a mere scale estimate, we not only achieved (confirmed) a simplification, but also 
conclude that the boundary layer thickness increases with a; 1 / 4 . 



2.8.2 Approximate equations 

We summarize the the approximate equations stating continuity, 

u x + w z = 

incompressibility, 

up x + wp z = 



(2.8.25) 
(2.8.26) 
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horizontal momentum balance, 

= -p x + fiu zz (2.8.27) 

and vertical momentum balance 

= -p g - pg- (2.8.28) 
Let the stream function ip be defined such that 

u = fa, w = -fa (2.8.29) 

Equation (2.8.26) means that in the direction of the local velocity, the density is constant. 
Thus the the density can only be a function of the stream function, i.e., 

p = p(fa (2.8.30) 

Let us check that (2.8.30) implies (2.8.26). Clearly 

dp dp 

px = -t-Vx, Pz = -TTVz, 
dip dip 

so that 

dp _ p x _ p y 

dip fa ip z 

Eq. (2.8.26) follows by using (2.8.29). 

Eliminating p between Eqn. (2.8.27) and Eqn. (2.8.28) we get 

pu zzz + gp x = 0, 

which can be rewritten in terms of the ip, 

vfazzz + -^fa = 0. (2.8.31) 
p dip 

Clearly this is a boundary layer equation. For a weakly stratifed fluid, we approximate p by 
a constant in the boundary layer. 
The boundary conditions are: 

u, ±oo. (2.8.32) 

which implies that 

u z , u zz , u ZZZ) ... > 0, z > ±oo 
It follows from Eqn. (2.8.31) that 

fa = w _> z -> ±oo (2.8.33) 

if dp /dip 7^ as z — > ±oo. Integrating Eqn. (2.8.25) from z = — oo to z = oo, and taking 
into account Eqn. (2.8.33), we get 

d 

— I u dz = 0. 

ox 
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Therefore, 

/oo 
udz = —Q (const.) (2.8.34) 
-oo 

where Q denotes the steady discharge. 

In terms of ip, the boundary conditions Eqns. (2.8.32), (2.8.33) and the integral constraint 
Eqn. (2.8.34) can be expressed as 

ipx, fa, tpzz -> (z -> ±oo) (2.8.35) 

and 

ip(x,oo)-ip(x,-oo) = -Q. (2.8.36) 

Remark: Unlike the 2-D laminar jet problem discussed before, the momentum flux rate 
is relatively unimportant in the present case of slow flow. 

2.8.3 Similarity solution for linear stratification 

In the present case u < hence ip increases as z decreases. Therefore dp /dip and dp/dz have 
opposite signs. For a stably stratified fluid, p increases as z decreases. Hence dp /dip > 0. In 
the special case where 

g dp 

— — = const. = v c > (J 
p dip 

(2.8.31) becomes linear 

ipzzzz + cip x = 0. (2.8.37) 
Let us look for a one-parameter transformation 

x = \ a x*,z = \ (3 z*,ip = X^ifj* 

such that Eqns. (2.8.35) - (2.8.37) are invariant. 

Since Eqn. (2.8.35) is homogeneous, it gives no information about the choice of a,/3, 7. 
From Eqns. (2.8.36) and (2.8.37), invariance of the boundary value problem requires that 

a = 4(3, 7 = 0. 

Therefore, a similarity solution can be found in the form 

ip £l s fc\ 1/4 
Q 

From Eqn. (2.8.38), it follows that 



/cX 1 / 4 

= f(v), V={-) z ( 2 - 8 - 38 ) 



w = -iP x = Q lf{rj) (2.8.39) 
u = ip z =^y 4 Qf'(ri), (2.8.40) 
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^ zz =^) 1/2 Qf"(v) (2.8.41) 

1>zzzz = -Qf""(v) (2-8.42) 

Ob 

Eq(2.8.37) then becomes 

4f""- v f' = (2.8.43) 
and Eqns. (2.8.35) and (2.8.36) become 

Vf',f,f"^0 fo-±oo) (2.8.44) 

/(oo) - f(-oo) = -1. (2.8.45) 

Remarks: 

(i) Since m(x, 0) must be finite, /'(0) is finite. It follows that 



v>(x,0)=-%-f'(v) 



= 0. 

77=0 



suggesting that w is odd in rj and antisymmetric about z = 0, in the x, z plane. Hence f{rj) 
is also anti-symmetric in z about z — 0, i.e., f{—rj) = —f{rj). Eqn. (2.8.45) may be repleced 
by 

/(oo) = -1/2 (2.8.46) 

It is always possible to take /(0) = 0, amounting to designating the axis as the streamline 
ip = 0. 

The anti-symmetry of w also implies symmetry of u in z about z — 0. These properties 
are the consequence of the assumption that dp /dip is constant and cannot be expected in 
general because of gravity. 

(ii) The similarity transformation can also be derived heuristically by examining Eqns. 
(2.8.35) - (2.8.37). Specifically, Eqn. (2.8.37) implies that 

-i- ~ — , thus 5 ~ x 1 / 4 
<5 4 x 

where 5 is the boundary layer thickness. Moreover, Eqn. (2.8.34) gives u ■ 8 ~ 0(1). 
Therefore, 

u ~ - ~ x 1//4 , and ^ ~ u5 ~ x°. 
o 

These results agree with (2.8.40). 
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2.8.4 Analytical Solution 

Let 

9(V) = f'(v)- (2-8.47) 

It follows from Eqn. (2.8.43) that 

Ag"'-r]g = 0. (2.8.48) 

Furthermore, in view of Eqn. (2.8.44), g and its derivatives vanish at rj = ±oo. Let us apply 
Fourier transform defined by 

1 r°° 

g{k) = — e- ik *>g(r))dr); 

Z7T J-oo 

Eqn. (2.8.48) is then converted to 

3~ , d 9 



The solution is 



Ait'g-i — 



g(k) = Ae 



where A is a coefficient. 

Taking the inverse Fourier transform of g, we get 



/OO P OO 

e -k*+ik v dk = 2 A cos k V e" fc4 dk. (2.8.49) 
-oo JO 

It is readily seen from Eqn. (2.8.49) that g{—rj) = g{i]), confirming the earlier arguments on 
the symmetry of u and f'(i]). 

In view of Eqns. (2.8.46) and (2.8.47), the constant A must be chosen such that 



/"CO POO 

2A drj dk cos kr) e 
Jo Jo 



1 

~ ~2 

It is shown in the Appendix that the double integral is — n/2 so that 



In particular, 

o(0) = 2A / e~ k dk = — ^— = -0.2885. 
v ; Jo 4r(3/4) 

The integral in Eqn. (2.8.49) may be evaluated numerically. Alternatively, Eqn. (2.8.48) 
is solved numerically by Chebyshev-polynomial method for unbounded domains (as, for 
example, J.P. Boyd (1992), J. Comp. Phys., 45, 43-79). Homework by Prof C. S. Wu, 
1.63 class of 1995, now at U of Wisconsin, Madison) 

From Figure 2.8.4 ( Prof. T.S. Yang, National Cheng Kung University, Taiwan (1.63 
Class of 1995) it is seen that, away from the center line of the sink, there are regions in 
which u > 0. Under the assumption that c > 0, we have, dp/dy > so that locally the fluid 
may be statically unstable. 



s 




Velocity profile of a 2D slow flow into a line sink in a density-stratified fluid. Calculations 
by G. D. Lee, 2002 

Appendix: A double integral (Yang, 1995, Homework) 

Consider the double integral 

I(n)= dke~ k l drj cos krj (2.8.50) 
Jo Jo 

Note that _ 

/oo 
5(k)e lk7] dk = 1 
-oo 

By inverse transform we get the Fourier integral representation of the delta function: 

5(k) = — / dke~ lk71 = - dr] cos krj 

Now we use this result in (2.8.50) to get 

I(n) = 7T / 5(k)e- k "dk = - (2.8.51) 
Jo 2 

It is interesting that the result is independent of n. Let us verify (2.8.51 ) by independent 
calculations for n — 1, 2: 



Since 

/■oo I 

/ dk e~ k cos krj = 

V2 



roo poo 

1(1) = / dke~ k / dr] cos krj 
Jo Jo 

/ dk e~ k cos kri = 

Jo 1 + ' 



and 



f°° 1 _-, oo 7T 

/ an = tan ri = — 

Jo 'l + m 'o 2 



(2.8.51) is proven. 

1(2) = / dke~ k / dr] cos 

JO JO 

Since 



and 




(2.8.51) is proven again. 



